Integrable dynamics of Toda-type 
on the square and triangular lattices 

P.M. Santini*, A. DoliwaJ and M. Nieszporski^ 
February 2, 2008 

Abstract 

In a recent paper we constructed an integrable generalization of 
the Toda law on the square lattice. In this paper we construct other 
examples of integrable dynamics of Toda-type on the square lattice, 
as well as on the triangular lattice, as nonlinear symmetries of the dis- 
crete Laplace equations on the square and triangular lattices. We also 
construct the r - function formulations and the Darboux-Backlund 
transformations of these novel dynamics. 

1 Introduction 

The Toda lattice [1, 2, 3] 
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where Amfm — fm+i ~ fm is the difference operator and qmit) is a dynamical 
function on a one dimensional lattice, is one of the most famous integrable 
nonlinear lattice equations. It describes the dynamics of a one - dimensional 
physical lattice, the masses of which are subjected to an interaction potential 
of exponential type. The infinite, finite and periodic Toda lattice (1), as well 
as its numerous extensions [4, 5, 6, 7, 8, 9, 10, 11, 12], have applications in 
various other physical and mathematical contexts [13, 14, 15, 16, 17, 18, 19]. 

Motivated by these results we find it important to construct integrable 
generalizations of the Toda law (1) to regular planar lattices; i.e., to the 
square, triangular and honeycomb lattices. Since the spectral problem as- 
sociated with (1) is an "integrable" discretization of the one-dimensional 
Schrodinger operator [20, 21] (where by integrable we mean that the op- 
erator admits, as its continuous counterpart, a large set of continuous and 
discrete symmetries, like the Laplace and Darboux transformations (DTs)), 
such a project requires the identification of proper integrable discretizations 
of self-adjoint second order operators on the plane first. A key progress in 
this direction was made in [22], where it was established that the self-adjoint 
scheme on the star of the triangular lattice admits Laplace transformations, 
and in [23, 24], where it was established that the self-adjoint schemes on 
the stars of the square, triangular and honeycomb lattices admit DTs, as 
their natural continuous counterparts. In addition, in [24], a novel discrete 
time dynamics on the triangular lattice was introduced, in connection with 
its Laplace transformation. To construct integrable nonlinear dynamics as- 
sociated with these self-adjoint operators, gauge equivalent to the discrete 
Laplace equations on weighted graphs, is the main goal of the paper. 

It is necessary to mention that these three planar schemes (on the square, 
triangular and honeycomb lattices) are directly connected (see [25] and [24]), 
via the sublattice approach [25] , to the so-called discrete Moutard [27, 29] (or 
S-quadrilateral [26]) lattice in Z^, and therefore they are all reductions of 
the multidimensional (planar) quadrilateral lattice [28, 30, 31, 32]. We also 
remark that the above three linear schemes are distinguished examples of 
Laplace equations on graphs, obtainable from the discrete Moutard equations 
on bipartite planar quad-graphs [33, 34, 35, 36]. 

Using the self-adjoint scheme on the star of the square lattice as spec- 
tral operator, we have recently constructed in [37] the following example of 
integrable Toda-type dynamics on the square lattice (together with its asso- 
ciated T - function formulation, its Darboux and Darboux-Backlund trans- 
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formations and some examples of explicit solutions): 

€m+l,7l+l 

(2) 

where qm,n{t) , im,n{t) are dynamical functions on the square lattice. 

Motivated by the above results, in this paper we construct and study 
other examples of integrable dynamics of Toda-type on the square lattice, 
as well as on the triangular lattice. In addition, we present their r function 
formulations, in which the r function of the BKP hierarchy [38] plays a central 
role, due to the already mentioned common origin of these schemes. The 
integrability of the dynamics in question manifests here in the construction 
of the Lax pair and Darboux - Backlund transformations (DBTs) . 

We remark that, due to the intimate connections between the self-adjoint 
schemes on the triangular and honeycomb lattices [24] , it is possible, in prin- 
ciple, to construct integrable Toda-type dynamics on the honeycomb lattice 
from those on the triangular lattice. This project will be developed elsewhere. 
Another interesting problem for future research is to establish connections 
between these Toda-like systems and the corresponding Lotka-Volterra sys- 
tems (see e.g. [39, 40]), as well as the connection, via the sublattice approach, 
between these Toda-like systems and the integrable dynamics on the discrete 
Moutard lattice introduced in [41]. 

The paper is organized as follows. In §2 we construct an integrable dy- 
namics of Toda - type on the square lattice, invariant under 7i/2 - rotation, 
its r - function formulation, and its two natural reductions transforming into 
each other under a 7r/2 - rotation. One of these two reductions coincides 
with the 2D Toda system (2) introduced in [37]. In §3 we construct an in- 
tegrable dynamics of Toda type on the triangular lattice, invariant under a 
n/3 - rotation, its r - function formulation and its natural reductions. The 
DBTs for all the above systems are presented in §4. 

2 Dynamics on the square lattice 

In this section we construct examples of integrable dynamics of Toda - type 
on the square lattice. To simplify the form of the equations,from now on, 
we will be using the following notation: / instead of /„^„, f±i instead of 
/m±i,n, /±2 instead of fm,n±u /±i±2 instead of fm±i,n±u f±i±i instead of 
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Figure 1: The square lattice and the points involved in the commutation. 

/m±2,n and /±2±2 instead of fm,n±2- Moreover we denote by Ti and T2 the 
basic translation operators acting on the lattice, i.e. Tj/ = f\, i = 1, 2. 

The wanted dynamics are associated with the linear self-adjoint 5-point 
scheme 

+ + + -B-2*-2 = (3) 

on the star of the square lattice, involving its black centre • and the four 
vertices of the star, denoted by the symbol □ in Fig 1. In equation (3) the 
eigenfunction is defined at the vertices of the graph, while the fields A, B 
are defined on the non oriented edges of the lattice. Equation (3), a natural 
discretization of the self-adjoint second order equation 

(a^^)^ + {h%)y = (4) 

admits, like its continuous counterpart, DTs [23]. 

We restrict our investigation to evolution equations for ^ involving only 
the 4 vertices □ of the 5-point scheme: 

— = CK^l + + 7^2 + (^*-2, (5) 

where the fields a,/?, 7,5, defined on the oriented edges of the lattice, will 
be specified in the following. A term proportional to ^ in (5) can always be 
expressed, using (3), in terms of the values of ^ at the 4 vertices □ of the 
star; therefore it is omitted. 
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We remark that, due to the 7r/2 - rotation symmetry of the square lattice, 
under which the two basic translations Ti,T2 transform as follows: 

T,^f,^ T2, T2^f2^ Tr\ (6) 

the coefficients of the 5-point scheme and of the evolution equation (5) are 
subjected to the following transformations: 

A^A^B, B^B^A_,, F^F^F, ^ 
a— >q; = 7, P ^ P — 6, 7— >7 = /3, S ^ 6 — a. 

The compatibility between equations (3) and (5) leads to an equation 
involving the values of ^ at all the marked points o, •, □, X, in Fig.l. 
Using the scheme (3) centered at the origin and at the points □, one ex- 
presses the values of ^ at the origin and at the points X in terms of the 8 
independent vahics of ^ at the points □ and o. As a result of this procedure, 
the compatibility condition becomes a linear equation for the 8 independent 
values of ^ at the points □ and o. Equating to zero their 8 coefficients, one 
obtains a determined system of 8 nonlinear equations for the 8 coefficients 
A,B,C,F,a,fl^,6. 

In the rest of this section we report the results of the analysis of such 
system leading to Toda-type dynamics. 



2.1 Rotationally invariant dynamics 

Setting 

a^^A, /? - 7 = -5, 5= --B_2, (8) 

2 2 2 2 ^ ^ 

where ^, rj are lattice fields to be specified, the corresponding evolution for ^ 

^ = 1 {A^l - + I (5^2 - 5-2^-2) (9) 

is compatible with the 5-point scheme (3) iff (up to a trivial gauge transfor- 
mation) the coefficients A, B, F, ^, 7] satisfy the following determined system 
of 5 nonlinear equations: 

f + 6^' - ^i^'-i + V2B' - V-2B'_2 = 0, 

if + lA,(eF)=0, + lA,(ryF) = 0, (10) 

AB{^ + ri)^ A2B^{^ + ri)u, A^Bi^ - 7^)2 = AB^{^ - ri),. 
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Equations (10b), (10c) suggest the introduction of the new fields q,r de- 
fined by 

dq dr ^ , , 

With this choice: 

A = ae^^i^ B = be'^^'^ (12) 

where a, b are arbitrary constants. Choosing, w.l.g., a,b— 1, the system (10) 
takes the form of the following novel integrable generahzation of the Toda 
law to the square lattice: 

(i+v) _ _AiA2(4^) _ AiA2(2^) 



Remark 1. In the natural 1-dimcnsional limit in which all the fields are 
invariant under the T2 translation, equations (13c, d) imply that ^ and r] are 
constant, and equation (13a) reduces to the 1-dimensional Toda lattice (1). 

Remark 2 Using (12) with a = b = 1, the 5-point scheme (3) takes the 
following form: 

Ti 1^ r 7^^r2^^^ r m/, ^^^^ 

r = e~2 , r = e"^. 

It is easy to verify that the spectral problem (14) reduces, in the natural 
continuous limit, to the stationary Schrodinger equation in the plane: "ifxx + 
+ ~ l^ therefore a natural integrable discretization of the 
Schrodinger operator, more general than that introduced in [23] . 

Remcirk 3 Using (6), (7), (8) and (11), it is easy to verify that, under a 7r/2 
- rotation, 

= V, V^V' = -^, q^q' = r, r ^ r' = -q; (15) 

from which it follows that the system (10) (or (13)) is invariant under this 
t r ansf ormation . 



6 



2.2 Reductions not invariant under rotation 

The system (10) (or (13)) admits two distinguished reductions for ^ = ±77. 
1. The reduction ^ = rj. In this case, the Lax pair (3), (9) reduces to 



(16) 



f = I (A*i - + S*2 - S-2*-2) . 

and the nonhnear dynamics (10) reduces to 

f + - ^-lAi, + - e-25^2 = 0, 

if + iA,(^F)=0, + iA2(eF)=0, (17) 

AB^ = ^2^162; 

Integrating equations (17b), (17c) and using (11), which imphes that r = 
g (r = r), one recovers the Toda type system (2), rewritten here, for com- 
pleteness, in the new notation: 

el (tl) - A, (ee-ie^--0 + A2 (ee-2e^--) , 

' 



and the associated 5-point scheme is the discrete Schrodinger equation 

r o _ r F _ p _ _2 (19) 



A=i^, B = f^, r = e-§, F- ^, 



introduced in [23]. 

2. The reduction ^ = —r]. In this case, the time evolution of ^ reads 

^ = I (A*l - - S*2 + 5-2*-2) , (20) 

at 2 

and the nonlinear dynamics (10) reduces to: 

f + ^^A^ - e-i^'_i - ^2B' + i-2B\ = 0, 

if + iAi(eF) = 0, -iA2(eF) = 0, (21) 



Equivalently, using (11) and noting that, in this case, r — —q (T — l/F), one 
obtains the Toda - type system: 



el (II) = Ai (a-ie^--0 - A2 (a-2e-^--) , 



(22) 



whose 5-point scheme is another variant of the discrete Schrodinger equation: 



(23) 



Wc end this section remarking that, due to the transformations (7), (15), 
the reduced systems (18) and (22) transform into each other under a 7r/2 - 
rotation. 

2.3 r - function formulations 

Motivated by the sublattice approach [25] for the self-adjoint 5-point scheme 
(3), we introduce two potentials r and f via equations 



B 



J2lL 
ff_i 



(24) 



These allow to resolve the algebraic part (10c) of the dynamic equations (10), 
with the fields ^ and r] expressed as follows 



2t2 



Tj 



T-lT-2 - ~-l-2T 

2t2 



(25) 



Then the remaining equations (10a) and (10b) form a system of three equa- 
tions for three fields r, f and F 



T_ir_2 r_i_2r 



0, 
0, 



2 dF 1 

+ 



n-2 ^ Ti 

dt ' ff_2 V f_2 T 
1 / ^-12 _ ^ 
ff_i V T^-l 



1 



T_iT_i_2 V T_i T_i_2 
1 / 7"_2-2 T_i_i_2 



T-2T_i_2 \ T_2 



T-1-2 
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Introduction of the fields q and r (or F and f ), which allowed to simplify 
equations (10b), suggests the introduction of yet other potentials h and h 
such that 



2 hih2 ^ I h 




It follows that 



^2 /il2^2 -A2 ^12^1 



p2 ^ ---^ p2 ^ (27) 

hhi hh2 

and that equations (10c) are identically satisfied, with the fields ^ and rj given 
as follows 

11 hhi2 hhi2 \ 1 / hhi2 hhi2 \ , . 

Moreover, equations (10b) reduce to two equivalent expressions for F 

hih2 d ( hi2\ hih2 d ( /i2 \ F 

log ^ = Tl— 17 1- = 17' (29) 



hhi2 dt\ h J hhi2 dt \ hi J 2 
while equation (10c) reads 

^ <j / hih2 d A A A ^ I ^^1 ^ ^2^12 I hih2 ( hjl ^ /tl2^2 j ^2^1 ^ 

dt\hhi2dt\ hJJ yhih hi2h2j^hhi2 \hhi h2hi2j_^hi2h 



hh2 hihi2\ /i2^i ( hi2hi h2h\ hih2 
h2h hi2hi J ^ hhi2 \h1h12 hh2 J _^ hi2h 



0. (30) 



Therefore the introduction of the potentials h and h allows one to rewrite the 
Toda-like system (10) as a coupled nonlinear system of two equations (the 
first equation of (29) and equation (30)). 



3 Dynamics on the triangular lattice 

In this section we construct some examples of integrable dynamics of Toda 
type on the regular triangular lattice. We recall that, on the triangular 
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Figure 2: The triangular lattice and the points involved in the commutation 

lattice, the three main translations Ti,T2,T3 in the directions 1, 2 and 3 are 
not independent, being connected by the relation 

= T2, (31) 

and hence /g = /_i2, f-3 = fi-2- 

The integrable dynamics of Toda type are associated with the linear and 
self-adjoint 7-point scheme 

A^l + + + 5-2*-2 + C*3 + C_3*-3 = F^, (32) 

on the star of the triangular lattice, involving the black centre • and the 
6 vertices denoted by the symbol □ in Fig. 2. In equation (32) the eigen- 
function ^ is defined at the vertices of the graph and the fields A, B,C are 
defined on the non oriented edges of the lattice. Equation (32), a natural 
discretization of the most general self-adjoint second order equation on the 
plane: 

(aV'x)x + {Hy)y + {ci^x)y + {c^y)x = /V', (33) 

admits, like its continuous counterpart, DTs [23]. 
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As in the previous section, we restrict our investigation to evolution equa- 
tions for ^ involving only the 6 points □ of the 7-point scheme: 

_ = a'^i + + 7^'2 + 5^-2 + e^'a + t^^'-s- (34) 
dt 

We remark that the regular triangular lattice possesses a 7r/3 - rotation 
symmetry, under which the 3 basic translations transform as follows: 

Ti^Ti=T2, T2^T2 = T3, T3^T3 = Tf\ (35) 

Correspondingly, the coefficients of the 7 - point scheme (32) and of the 
evolution equation (34) transform as follows: 

A^B, B^C, C ^ 

a— >q; = 7, P ^ P — 5, 7— >7 = e, 5^6 — 1/, (36) 
e— >e = /3, v ^ i> — a. 

We proceed adopting the same strategy as in the previous section. The 
compatibility between equations (32) and (34) leads to an equation involving 
the values of at all the 19 marked points in Fig. 2. Using the scheme (32) 
centered at the origin and at the points □, one expresses the values of ^ at 
the origin and at the points X in terms of the 12 values of "if at the points 
□ and o. As a result of this procedure, the compatibility condition becomes 
a linear equation for the 12 independent values of ^ at the points □ and o. 
Equating to zero their 12 coefficients, one obtains an ovcrdetermined system 
of 12 nonlinear equations for the 10 coefficients A, B,C, F,a, l3,'~f,6,e,i'. It 
turns out that, due to the relation (31) among the three main shifts, such 
overdeterminacy is resolved, and one construct integrable nontrivial dynam- 
ics. 

In the rest of this section we report the results of such analysis, leading 
to the Toda-type dynamics on the triangular lattice. 

3.1 Rotationally invariant dynamics 

Setting 

a^^A, 13= 7 = -5, 5= --B_2, e = ^C, u= --C_3, (37) 

2 2 2 2 2 2 ^ ^ 
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where ^, 77, C a-re lattice fields to be specified, the corresponding evolution for 

^ = I (A^l - + I (5^2 - 5_2^_2) + I (C^3 - C_3^-3) (38) 

at Z Z Z 

is compatible with the 7-point scheme (32) iff the coefficients A, B, C, F, ^, 77, ^ 
satisfy the following determined system of 7 nonlinear equations: 

f + iAi {ii-,A\) + iA2 (w_25!2) + ^As (CC-aC^s) = 0, 
if + \^^m + O2 + l^iv + 0-3 = 0, 

if + l^^ivF) + i^(e - C)3 - |^(e - C)i = 0, 

if + iA3(CF) + i^(e + r/)-i - i^(e + r7)2 = 0, (39) 
ASi(e- 77)1 = ^2^(^-77)2, 

AC7(C + C)=A3Ci(e + C)2, 
S_3Ci(77-C)l = SC_3(77-C). 

We remark that, due to equations (39e)-(39g), equations (39b)- (39d) can be 
rewritten in the following conservation-like form: 

|(ln A^) + A,m - A2 (^(e - V)) - A3 (^(e + 0-3) = 0, 
l(lnS^) + A2(77F) + A3 (^(77 - 0) - Ai (^(e - 77)) = 0, 

l(lnC^) + A3(CF) - Ai (^(e + C)-i) - A2 {^{V - 0) = 0. 

(40) 

Remcirk 4 Under the transformation (36), the coefficients ^,77,C transform 
as follows: 

= v^v^c C^C = -e, (41) 

and. as it is easy to verify, the nonlinear system (39) is invariant under a 7r/3 
- rotation. 
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3.2 Reductions not invariant under rotation 



The nonlinear system (39) admits the reductions = [rj = Q and [C = 
—^], and the following combinations of them: = rj, rj = = rj, ( = 

— and [rj = (, ( = —C,]. They give rise to six integrable dynamics on 
the triangular lattice. It follows that these dynamics are not rotationally 
invariant, but they transform one into the other in the way summarized in 
Fig. 3. 




Figure 3: The bold arrows describe the 6 reductions of the rotationally in- 
variant system. The dashed arrows describe how such reductions transform 
one into the other, under a 7r/3 rotation. 

We write down explicitely the two reductions = rj] and = r], r] = (], 
since all the others can be generated from them through rotations. 
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The reduction i — rj. In this case, the evolution of ^ reads: 

at I I 
and the nonhnear system (39) reduces to the six equations 

f + |(Ai {iUA\) + {iUBl,)) + iAa {CUCl,) = 0, 



d(lnA^) 
dt 



+ Ai(eF) - A3 + 0-3) = 0, 

^ + A,(eF) + A3(^(e-C))=0, 

^ + A3(CF) - A, (^(e + 0-1) - (^(e - 0) = 
^C(e + C) = ^3Ci(e + C)2, 

s_3Ci(e-C)i = 5C'_3(^-C)- 

The reduction {i = r], rj = (). In this case, the evolution of ^ reads: 

^ = I {A^i - + - S_2*-2 + C*3 - C_3*-3) , 

and the nonlinear system (39) reduces to the five equations 

+ Ai (a-l^'-l) + A2 (a-25^2) + As (CC-3^^^3) = 0, 

^ + Ai(eF) - 2A3 (^-^-^^-3) = 0, 

« + A2(eF)=0, 

^ + A3(eF) - 2Ai (^e-i) = 0, 
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In analogy with the previous examples on the square lattice, equations 
(45b) - (45d) suggest the introduction of the new fields g, p, a defined by 



With this choice: 

where a, 6, c are arbitrary constants. Choosing, w.l.g., a,b,c = 1, the system 
(45) can be rewritten in the following Toda-like form: 

^1 (11) + ^1 (ee-ie^i«-+^--'--) + A2 (ee-2e^^''-) + 

dE ^ _2^gA3(«+f)-iAia_i^ (48) 



da 
dt 



-2^e^i(«+t)-5^3''-3, 



^2 = ^e<13-1'2+<ll-<l+^i^'i(^-l+^3P-3) _ 



3.3 Reductions to dynamics on the 7? graph 

We remark that the reduction ^ = t], 7] = ( is compatible with the condition 
B = 0, for which all the connections in the direction 2 are broken, and the 
triangular lattice reduces to the rhombic lattice in Fig.4a. Then the direction 
2 should be renamed 13 and the system (45) becomes the integrable system 
(17) on the rhombic lattice of Fig. 4a (on the graph). 

Analogously, it would be possible to show, f.e., that the reduction ^ = 
—Tj = — C is compatible with the condition C = 0, for which all the connec- 
tions in the direction 3 are broken, and the triangular lattice reduces to the 
rhombic lattice in Fig.4b. The system obtained in this case is the integrable 
system (21) on such a rhombic lattice. 
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Figure 4: a) The rhombic lattice for B = 0. b) The rhombic lattice for 
C = 0. 



3.4 r - function formulations 

Motivated by the sublattice approach [24] for the self-adjoint 7-point scheme 
(32), we introduce three potentials r, f and f such that 

A=^^, B = ^, C=^. (49) 

Tf2 Tsfs 

The algebraic part (39e)-(39g) of the nonlinear system (39) is then resolved 
by the parametrization 

2r^ 
1 

V = TT^Crsr -f2T_i + ff3), (50) 

( = ^(^3^ - ^2^-1 - rfg). 

The remaining part of the system gives a system of four equations for the 
four fields F, r, f and f. 



4 Darboux - Backlund transformations 

In this section we present the DBTs for the above Toda-type systems. 
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4.1 DBTs for the Toda-type system on the square lat- 
tice 

The Lax pair (3), (9) is covariant under the Darboux transformation 

{^,A,...,r))^{^,A,...,'n) (51) 
given by the hnear system 

A2 (k^) = -A_,e9.,A_, (I) , 
Ai [K^) = 5_2^^-2A_2 (f ) , 



^-1^-2^-1-8-2 



(e)-2 ie)-i 



+ 



(52) 



^^_l_2A-lS-i_2(i^ [(|)-l-2-| 

where ^ is a particular solution of (3), the gauge function K must obey 

1 '^^ + (^-1^-1 + ^-2^-2)^ - (^^-2^-1-2 + ^-li?-l-2)%?^ + 



K dt 

(^_l_2A_i_2 + ^S_2) 



2^ 461 y^-^^'—i--^ I '^-i-^-i-zy 46l_i_2 
^-(M_i + e-l-2S_,_2)(a 



0, 



(53) 



and the transformation of the other fields is given by 



B = 
F = K 

i+v 



KK2 



2/ 1 



+ 



+ 



^A_i5_i_2(^-r?)-i, 
^A_iS_2(e + r/). 



+ 



-1-2(7-1-20 



-I 



(54) 



So the formulae (54) are the Backlund transformations (BTs) for the Toda- 
type system (10) on the square lattice, i.e. A, B, F, ^, is a new solution of 
(10). On the level of the r-functions the transformation is given as follows 



T I— > tiT. 



(55) 



We remark that, for C, = ^i], the above transformations become the DBTs 
for the reduced systems (17) and (21). 

The spatial parts of the above DBTs were already written in [23]; the 
temporal parts, describing the time dependence of the transformed solution 
^, and the transformation law for the coefficient ^, 77, C, are new ingredients 
of this paper. 
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4.2 DBTs for the Toda type system on the triangular 
lattice 

First, for aesthetical reasons, we introduce function 

S := C_2. 

The Lax pair (32), (38) is covariant under the Darboux transformation 

(*,A---,C)^(^,A---,C) (56) 

given by the hnear system 
l(ir^) = -i{^_i_2A_i5_i_2(e - r/)_ivl>+ 

^_2A_i_i5(e + 0-1^-1-1 - e.,B_2.2S{v - 0-2^-2-2 + 
f-f;';' {V - 0-1-2 + ^-l-2^A-l-2(e - V)-2 + ^-2A-l-2i?-l-2(e + C)-l-2]^-l- 
+ 0-1-2 - ^-l-2-5S_i_2(e - + ^-lA_i_25_i_2(r/ - 0-l-2]^-2- 

[^%^{^ - VU - ^-25S_i_2(e + 0-1 + ^-l5A_i_2(77 - 0-2]*-l-2}, 

(57) 

where ^ is a particular solution of the Lax pair (32), (38), P is given by 

P := M_iS_2 + e^iA^iS + 
and K is given by the quadrature 

1 dK _ zlLS Pf) B~l-2 (c _ r ^-1 ^-1-1 B^2 _|_ B^2-2 S-i _|_ 5^2 1 

K dt 4 i-^'^-i-2 B-2 'l'-^ye-2P e-i-2P-i 9~iP ~'~ 9-1-2P-2 e-iP-1 e_2P-2J 

D /3 -4_i_2/j^ I A\ r -4-1 ^-1-1 I -B-2 -5-2-2 I S-l 5-2 1 I 

^-1^-2X1777 + Cj-i-2[el7P - e-i-2P-i + " 0-i_2P_2 + eHi^ ~ er;iriJ + 

P -B-1-2 / _ A\ r A-i _|_ P-2 B-2-2 5-1 _|_ 5-2 u 

-^-2'^-1B_2_2 V' ''/-1-2L6»_2P 6'_1_2P-1 0-lP 0_l_2P-2 e»-lP-l 0-2P-2JJ' 

(58) 

The new eigenfunction ^ is a solution of the Lax pair (32), (38) with the 
new coefficients 

A _ KK^A , 

^ ~ e-2P^-i' 

R _ KK2 p 

- e_i_2P_i_2'^-i-2' 
F = K'^ {j-^ + e_i_2P-i el^P 6>_i_2P-2 6>_iP_i + el^Pl^) (^^) 

e-r) = ^WiV'^ (e-r^)-i, 

^~ + C = %^(^ + 0-i, 
^~-C = %l|^(^-0-2- 
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Therefore formulae (59) constitute the BTs for the Toda type system (39) 
on the triangular lattice. 

We would like to mention that the first two equations of (57) can be easily 
inverted 

In addition, equation (58) and the fourth equation of (59) can be rewritten 
by means of " new" solutions as follows 



d f 1\ I ~ A A_i\ ( B 5_2 \ z S2 Si 



dtXKj 2 ^ \Ki K_i ' \K2 K_2 \ i^i- 




F 

: i.e., is eigenfunction of tilded Lax pair. 

On the level of the r-functions, the transformation is given as follows 

T AT_i, T y-2T-2, T ^ : • (olj 

V ^-l J -2 

As before, the DBTs (59) are consistent with all the reductions of the 
Toda type system (39). 

The spatial parts of the above DBTs were already written in [23]; the 
temporal parts, describing the time dependence of the transformed solution 
^, and the transformation law for the coefficient C,ri,(, are new ingredients 
of this paper. 
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